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Abstract 

A computational procedure is presented for the 
geometrically nonlinear analysis of aircraft tires. 
The Space Shuttle orbiter nose-gear tire was mod- 
eled through use of a two-dimensional laminated 
anisotropic shell theory with the effects of variation 
in material and geometric parameters included. The 
four key elements of the procedure are (1) semi- 
analytic finite elements in which the shell variables 
are represented by Fourier series in the circumferen- 
tial direction and piecewise polynomials in the merid- 
ional direction; (2) a mixed formulation with the 
fundamental unknowns consisting of strain parame- 
ters, stress-resultant parameters, and generalized dis- 
placements; (3) multilevel operator splitting to ef- 
fect successive simplifications and to uncouple the 
equations associated with different Fourier harmon- 
ics; and (4) multilevel iterative procedures and re- 
duction techniques to generate the response of the 
shell. Numerical results of the Space Shuttle orbiter 
nose-gear tire model are compared with experimen- 
tal measurements of the tire subjected to inflation 
loading. 

Introduction 

Because of the axial symmetry of undeformed 
tires, it is desirable in their modeling and analysis to 
exploit the substantial capability that currently ex- 
ists for the numerical analysis of shells of revolution. 
The most commonly used approach for the analysis of 
shells of revolution is based on the representat ion of 
the shell variables and loads by a Fourier series in the 
circumferential coordinate 6 , combined with the use 
of a numerical discretization technique (such as finite 
elements, finite differences, or numerical integration) 
in the meridional direction. (See, for example, refs. 1 
to 7.) Such an approach has the major advantages 
of accuracy and stability (i.e., no locking or spurious 
modes) over approaches which use two-dimensional 
shell elements. Moreover, for linear problems of shells 
with uniform circumferential properties, the Fourier 
series representation permits separation of variables 
and the equations uncouple in harmonics. However, 
when applied to the analysis of tires the approach has 
the following drawbacks, which can make the compu- 
tational cost of the geometrically nonlinear analysis 
of tires quite expensive: 

1. For geometrically nonlinear problems, the 
unknowns associated with different harmonics are 
coupled (see refs. 8 and 9). 

2. Even for linear problems, because of the 
anisotropy of the cord-tire composites, the symmetric 
and antisymmetric responses (with respect to 0 = 0) 
asssociated with each harmonic are coupled. 


3. For localized loading (e.g., contact pressure on 
the tire), a large number of harmonics are needed to 
accurately predict the response. 

Research on tire modeling and analysis at NASA 
Langley Research Center has focused on developing 
accurate and cost-effective strategies for predicting 
tire response. Included in the research is the develop- 
ment of analysis procedures for substantial reduction 
of the computational expense resulting from the har- 
monic and anisotropic couplings (items 1 and 2) and 
the generation of the response associated with large 
numbers of harmonics (item 3). The present paper 
summarizes the status of these development activi- 
ties. To demonstrate the capabilities of the analysis 
techniques, numerical studies were conducted with 
the Space Shuttle orbiter nose-gear tire. Numerical 
results are presented for an inflated Space Shuttle 
tire, and these results are compared with experimen- 
tal measurements. 

Notation 

Aj j , B 2 j , Djj stiffness coefficients of two- 

dimensional shell model of tire 

(ij = 1,2,6) 

01 , 02 parameters defining geometry of 

tire cross section (see fig. 3) 

c t j stiffness coefficients of indi- 

vidual unidirectional plies 
(i, j = 1,2, 4, 5, and 6) 

d diameter of tire nylon cord 

{ E } vector of strain parameters for 

shell model of tire (sec eq. (2)) 

E Cl E r Young's moduli for nylon cord 

arid rubber 

{E} n vector of strain parameters 

associated with nth Fourier 
harmonic for shell model of tire 
(see eq. (2) and table 1) 

Ej o elastic modulus (see tables 7 

and 8 and fig. 9) 

E\, E 2 elastic moduli in direction of tire 

cord and normal to it 

epi cord end count, ends per inch 

fc volume fraction of nylon cord in 

individual plies of tire 

{/}■"■ vector defined in equations (5) 

G c ,G r shear moduli for tire cord and 

rubber 



{G}( n ) vector of nonlinear terms associ- 

ated with nth Fourier harmonic 
(see eqs. (5)) 

G 12, G 13 shear moduli in plane of tire cord 

G23 shear modulus normal to plane of 

tire cord 

{H} vector of stress-resultant param- 

eters for shell model of tire (see 
eq. (3)) 

{H} n vectors of stress-resultant param- 

eters associated with nth Fourier 
harmonic for shell model of tire 
(see eq. (3) and table 1) 

h total thickness of tire 

h nondimensional thickness of tire 

(see fig. 2) 

hk thickness of individual layers of 

two-dimensional shell model 

h G total thickness of tire at £ = 0 

(see fig. 3) 

[K], [A"] matrices defined in equations (6) 

[A"](^) matrix containing nonlinear 

terms (see eqs. (9) and (11)) 

[K]( n } linear matrix associated with nth 

Fourier harmonic (see eqs. (5)) 

[K 0 ] y [ K a ] submatrices of [A"]^ (see 

eqs. (8)) 

[A;], [k\ matrices associated with reduced 

equations (see eqs. (16) and (17)) 

l number of segments along tire 

surfaces 

M s , Mq, M s q bending arid twisting stress 

resultants (see fig. 1 and table 1) 

m number of displacement nodes in 

element 

N number of Fourier harmonics 

which are greater than or equal 
to 1 

N l shape functions used in approxi- 

mating generalized displacements 
and external loading 

N l shape functions used in approx- 

imating stress resultants and 
strain components 


N,,N 0 ,N s6 


n 


n 


Tin 


{P)AP} 

{P}M 

Pjl 

Po 

Po 

PsiPOiP 

Qs, Qo 
{«},{?} 

R\, R2 


r Q 


extensional stress resultants (see 
fig. 1 and table 1) 

vector normal to reference surface 
of tire 

Fourier harmonic (circumferential 
wave number) 

Fourier harmonic at which 
global approximation vectors are 
generated 

right-side vectors in 
equations (13) 

consistent load vector associated 
with nth Fourier harmonic for 
tire (see eqs. (5)) 

normal pressure components 
associated with nth Fourier 
harmonic for shell model of tire 
(see eq. (1)) 

intensity of localized normal 
loading of tire (see fig. 10) 

inflation pressure acting normal 
to inner surface of tire (see figs. 3 
and 4) 

intensity of external loading in 
coordinate directions (see fig. 1) 

transverse-shear stress resultants 
(see fig. 1 and table 1) 

load vectors associated with 
reduced equations (see eqs. (18) 
and (19)) 

principal radii of curvature in 
meridional and circumferential 
directions of reference surface of 
shell model 

normal distance from tire axis to 
reference surface (see fig. 3) 

= r at £ = 0 

submatrices of [A']( n ) (see 
eqs. (8)) 

meridional coordinate of tire (see 
fig. 1); number of parameters 
used in approximating each stress 
resultant and strain components 
in element 

polygonal arc lengths of segment 
of tire surfaces (see eq. (A2)) 


2 


u 

u, v, w 


W 


X 


X, z 


^3 


y; 


{Z}n 

a 

[r] 

712 ? 713 » 723 


E s ,£q, 2e s Q 


2£s3? 2^3 


£\ > £2 


9 

6 


O k 


K Si 2 k s Q 
«li«2 


strain energy density of shell 
model of tire 

displacement components of ref- 
erence surface of tire in merid- 
ional. circumferential, and normal 
directions (see fig. 1 and table 1) 

vector of generalized nodal 
displacement coefficients for tire 
(see eq. (4)) 

vector of generalized nodal dis- 
placement coefficients associated 
with nth Fourier harmonic for 
shell model of tire (see table 1) 

position vector for points to be 
interpolated along segment of tire 
surfaces (see eqs. (Al)) 

Cartesian coordinates (see fig. 3) 

coordinate normal to tire refer- 
ence surfaces (see fig. 1) 

position vector for data points of 
segment (see eqs. (Al)) 

vector of unknowns associated 
with nth Fourier harmonic (see 
eqs. (7)) 

contact angle in meridional 
direction (see fig. 10) 

matrix of global approximation 
vectors (see eqs. (14)) 

shear strain components in 
principal material axis 

extensional strains of reference 
surface of tire 

transverse-shear strains of shell 
model of tire 

in-plane strain components in tire 
cord direction and normal to it 

circumferential (hoop) coordinate 
of tire (see fig. 1) 

orientation angle used in equa- 
tion (24) and table 6 

orientation angle of tire cord, deg 

bending strains of tire 

principal curvatures in meridional 
and circumferential directions of 
reference surface of shell model 


« 1 . 0 ? ^ 2,0 
A 


Vc,V T 


v\2,v 2 \ 


o 


o 




principal curvatures at £ = 0 

tracing parameter identifying 
coupling between different 
Fourier harmonics (see eqs. (12)) 

Poisson's ratios of tire cord and 
rubber 

major and secondary Poisson's 
ratios in individual plies 

dimensionless coordinate along 
meridian (see fig. 2) 

tension factor 

normalized tension factor (see 
eq. (A4)) 

in-plane stress components in tire 
cord direction and normal to it 


ri2, Ti 3 , 7"23 shear stress components in 
principal material axes 


0 


angle between normal to merid- 
ian and x-axis (see fig. 3) 


< p$,(t>ft rotational components of refer- 

ence surface of tire (see fig. 1 and 
table 1) 

{'F} n vector of amplitudes of global 

approximation vectors (see 
eqs- (14)) 

Superscripts: 

i index of shape functions for 

approximating generalized 
displacements and external 
loadings; ranges from 1 to m 


l index of shape functions for 

approximating stress resultants 
and strain components; ranges 
from 1 to s 


r 

t 


Subscripts: 

max 

n 


number of iterational cycles 
matrix transposition 

coefficient of sine terms in 
Fourier series 

maximum value 
Fourier harmonic 


Analysis 


Mathematical Formulation 

In the present study we model the Space Shut- 
tle orbiter nose-gear tire using a moderate-rotation 
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Sanders-Budiansky shell theory with the effects of 
transverse shear deformation and laminated 
Anisotropic material response included (refs. 10 and 
11). A total Lagrangian formulation is used and the 
fundamental unknowns consist of the five generalized 
displacements, the eight stress resultants, and the 
corresponding eight strain components of the mid- 
dle surface. The sign convention for the different 
tire stress resultants and generalized displacements 
is shown in figure 1. The concepts presented in the 
succeeding sections can be extended to higher or- 
der shear deformation theories as well as to three- 
dimensional continuum theory. 

Spatial Discretization of the Tire 

Each of the generalized displacements, the stress 
resultants, and the strain components is expanded in 
a Fourier series of the circumferential coordinate 0. 
The discretization in the meridional direction is per- 
formed through the use of a three-field mixed finite- 
element model. The following expressions are used 
for approximating the external loading, strain com- 
ponents, stress resultants, and generalized displace- 
ments within each element: 


the quantities associated with the Fourier harmonic 
n. Note that the degree of the polynomial shape 

functions is lower than that of NK Moreover, 
the continuity of the strain components and stress 
resultants is not imposed at the interelernent bound- 
aries and, therefore, the strain and stress- result ant 
parameters can be eliminated on the element level. 

In equations (1) to (4), the range of i is 1 to rn, 
the number of displacement nodes in the element; 
the range of / is 1 to s, the number of parameters 
used in approximating each of the strain components 
and stress resultants. The shell variables without a 
bar are the coefficients of the cosine series, the shell 
variables with a bar are the coefficients of the sine 
series, and a repeated superscript denotes summation 
over its entire range. Henceforth, the vectors of 
the 10 generalized displacement parameters, of the 
16 strain parameters, and of the 16 stress-resultant 
parameters chosen in association with the harmonic n 
are denoted by {X} n , {// } n , and {E} n , respectively. 
These vectors can be decomposed into symmetric and 
antisymmetric sets (with respect to 9 = 0) as shown 
in table 1. 


p(s, 9) = N l I p l n cos n9 + pjj sin n0 J ( 1 ) 
Vn=0 n = 1 / 


/ UU CXJ \ 

E(s , 0) = V E '\ t cos n0 + En sin n0 J (2) 


v7l = 0 


n= 1 


Governing Equations 

The governing discrete equations of the tire are 
obtained through application of the three-field Hu- 
Washizu mixed variational principle (see ref. 12). 
If the number of terms (harmonics) retained in the 
Fourier series is N + 1, then the governing equations 
can be written in the following compact form: 


( OO DC 

Hi j cos ri0 -f h\ x sin 7i0 | (3) 

71=0 71 = 1 / 


X(S, 0) = N l I X n cosn0 + Y, X n Sin n6 ( 4 ) 

WO n=l / 


where N l are the shape functions used in approx- 
imating the generalized displacements and external 

loading in the meridional direction; N l are the shape 
functions used in approximating the strain compo- 
nents and stress resultants; p l n and p\ x refer to the 

normal pressure coefficients; X l n and lC n refer to the 
generalized displacement coefficients; E l n and E l n re- 
fer to the strain parameters; H l n and H l n refer to the 
stress- resultant parameters; and subscript n refers to 


' / (0) ’ 




' Zo ' 

/(i) 


kW 


Z\ 


> = 


< 

> 

. . 


K(X)_ 


- Zpj j 


( G^{Z Q ,Z X ,. 

..,Z N ) 1 


r p{ 0) ' 

G^(Z q ,Z u . 

■ ■,Z N ) 

I 

p( 1) 

+ < 


. - < 

> 

[g^'\z Q) z u . 

■ -,Z N ) , 

: 

. P( N ) ; 


where {Z) n (n = 0,1,.,., AT) is the vector of 
unknowns associated with the Tith harmonic and 
includes vectors of strain parameters { E } n , of stress- 
resultant parameters {H} n , and of generalized 
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displacements { X } n ; [K]^ is a linear matrix; {G}^ 
is a vector of nonlinear terms; and (P}( n ) is a con- 
sistent load vector. 

The following observations can be made about the 
governing equations (eqs. (5)): 

1. The first matrix on the left side of equations (5) 
is block diagonal, which is a direct consequence 
of the orthogonality of the trigonometric func- 
tions. The orthogonality of trigonometric func- 
tions leads to uncoupling of the equations associ- 
ated with the different Fourier harmonics for the 
linear case. For the nonlinear case, the vector 
{G}^ couples the unknowns associated with all 
the harmonics. (See, for example, ref. 9.) 

2. The contributions of the different Fourier harmon- 
ics and the anisotropic (nonorthotropic) mate- 
rial coefficients to the governing equations can be 
identified as follows: 


between the symmetric and antisymmetric 
shell parameters (see ref. 13 and table 2). 

3. If the vector {Z} n is partitioned into subvectors 
of parameters of strains, stress resultants, and 
generalized displacements, that is, 


{Z}n 




(7) 


then the matrix can be written in the 

following form: 


[K]W = 


A f > T A a 
-R f 


-R 

S l 0 + nS f 


So + nS 


(8) 


a. Fourier harmonics — The block-diagonal ma- 
trix [K]^ (n > 1) in equations (5) is linear 
in the Fourier harmonic n. Therefore, [K]^ 
can be expressed as the sum of two matrices 
as follows: 

[tf]W = [K] + n\K } (6) 

where both [K] and [K] arc independent of n. 
The nonlinear vector {G}( n ' is quadratic’ in n. 

b. Anisotropy (nonorthotropy)-- A unique fea- 
ture of the mixed formulation used herein is 
that the anisotropic (nonorthotropic) mate- 
rial coefficients are included only in the lin- 
ear matrix [AT]^. For the linear case, these 
anisotropic coefficients result in the coupling 


where the submatrices [K 0 ] and [K a ] contain the 
contributions of the orthotropic and anisotropic 
(nonorthotropic) material coefficients. The ex- 
plicit forms of the submatrices [K G \, [K a ], [5 0 ], [5], 
and [R] are given in references 14 and 15. 

4. The nonlinear vector {G}^ contains bilinear 
terms in {H} n and {A r } n as well as quadratic 
terms in { JVT } 7 , . 

Generation of the Nonlinear Response of 
the Tire 

For a given external loading, the governing non- 
linear equations (eqs. (5)) are solved by using the 
Newton-Raphson iterative technique. The recursion 
formulas for the rth iterational cycle are 



A' 11 ’ I 


'K m ~K {m) ... 

\ 

{r) 

' AZ„ ' 

(r) 



A' 1 ') 


o ... k" N) 



AZ, 


/ (1 > 



+ 



< 


► = - < 

> 


A'(*>_ 


0 . 

) 


. A Z.v . 


. / (lV) . 


and 


r ' 

(r+l) 

' Zo ' 

(r) 

r AZ 0 1 





AZ, 


► = < 


> + < 

> 



. Zn > 


k A Z N j 


( 10 ) 
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where 

\W J] = A] {G}{I) (I , J — l to N) ( 11 ) 

For each Newton-Raphson iteration (represented by eqs. (9) and (10)), another iteration loop is performed 
using the preconditioned conjugate gradient (PCG) technique to account for the coupling between the different 
harmonics (i.e., the submatrix i\ ' . i n the inner iteration loop the following uncoupled equations are 

solved: 

([A'] (0) + [i?] (0())(r) ) {A Z} ( Q r) = — {/} (0)(r) - A([F]( 01 ){AZ}i 

+ [¥](° 2 >{AZ 2 } + • • ■ + [F](° jV ){AZ N }) (r) 

[A'P^AZ}^ = - a([K]( 10 ){AZ} 0 + [K]W{AZ} 2 

+ - + 1 7f]< 1 '' tf| {AZ) w ) (r> 

= - {/|' ,V|,r| - X ([3?|< m >{ AZ} 0 + [3f](' vl *{AZ) 1 + . . j" 



where A is a tracing parameter which identifies the 
coupling between the different Fourier harmonics. 
When A = 1 equations (12) are equivalent to equa- 
tions (9), and when A = 0 the equations uncouple 
in harmonics. Note that because of the special struc- 
ture of the Jacobian matrix in equations (9), only the 
left side associated with the zeroth harmonic needs 
to be updated in each iteration. An efficient tech- 
nique is described in the next subsection for solving 
equations (12). 

Efficient Generation of the Response 

Associated With Different Harmonics 

An efficient procedure is presented herein for gen- 
erating the tire responses associated with different 
harmonics (solution of eqs. (12)). The basic idea of 
this procedure is to approximate the tire response 
associated with the range of Fourier harmonics, 1 < 
n < TV, by a linear combination of a few global ap- 
proximation vectors that are generated at a particu- 
lar value of the Fourier harmonic within that range. 
The full equations of the finite-element model are 
solved for only a single Fourier harmonic, and the re- 
sponses corresponding to the other Fourier harmon- 
ics are generated using a reduced system of equa- 
tions with considerably fewer degrees of freedom. 
The proposed procedure can be conveniently divided 
into two phases: (1) restructuring equations (12), for 
1 < n < TV, to delineate the dependence on the 
Fourier harmonic n, and (2) generating global ap- 
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proximation vectors (or modes) to approximate the 
response associated with a range of values of the 
Fourier harmonic and determining the amplitudes of 
the modes. Application of the procedure to stress 
and vibration problems of anisotropic shells of rev- 
olution is described in references 16 and 17. Its ap- 
plication to the solution of equations (12) is outlined 
subsequently. 

Restructuring of the governing equations . 
If equations (6) and (8) are used, the governing 
equations for the harmonic n (1 < n < N) can be 
embedded in a single-parameter family of equations 
and written in the following compact form: 

([K] + n[K}) {A Z} n = {P(n)} + A{P(n)} (13) 

The two vectors {P(n)} and {P{n)} are quadratic 
in n. 

Basis reduction and reduced system of 
equations . The basis reduction is achieved by ap- 
proximating the vector {A Z} n , for a certain range of 
Fourier harmonics, 1 < n < TV, by a linear combina- 
tion of a few global approximation vectors which are 
generated at a particular value of the Fourier har- 
monic within that range. The approximation is ex- 
pressed by the following transformation: 

{A Z}n = [r]{*} n ( 14 ) 


where [r] is a transformation matrix whose columns 
are the preselected approximation vectors, and {^} n 
is a vector of unknown parameters representing the 
amplitudes of the global approximation vectors for 
the harmonics n. The number of components of 
{4>} n is much less than the number of components 
of {A Z} n . 

A Bubnov-Galerkin technique is now used to re- 
place the original equations (eqs. (13)) by the follow- 


ing reduced equations in {'Fjn: 

({k]+n[k})m n = {q} + X{q} (15) 

where 

[k] = [if [£][ r] (i6) 

\k\ = [if[£][r] (17) 

M = [r ]‘{P(n)} (18) 

{?} = [lf{P(n)} (19) 


Selection and generation of global approxi- 
mation vectors . The global approximation vectors 
are selected to be the response associated with a sin- 
gle Fourier harmonic n (} and its various-order deriva- 
tives with respect to n. Henceforth, the derivatives 
of the response with respect to n are referred to as 
“path derivatives.” The matrix [T] in equations (14) 
is therefore given by 

[r] = [(AZ} l (AZ| £, {A Z) ...]_ 

(20) 

The path derivatives are obtained by successive 
differentiation of the governing equations (eqs. (13)). 
The recursion relations for the first three global 
approximation vectors can be written in the following 
form: 


Note that the left-side matrix in equations (21) to 
(23) is the same, and therefore it needs to be decom- 
posed only once in the process of generating all the 
global approximation vectors. 

Comments on proposed procedure . The fol- 
lowing comments are made concerning the forego- 
ing procedure for generating the responses associated 
with different harmonics: 

T The particular choice of the global approximation 
vectors used herein provides a direct quantitative 
measure of the sensitivity of the different response 
quantities of the tire to the circumferential wave 
number (the Fourier harmonic) n. 

2. For problems requiring large numbers of Fourier 
harmonics (e.g., 100 or more), the range of n is 
divided into intervals of fewer (e.g., 7) harmonics 
each; the global approximation vectors and re- 
duced equations are generated at an intermediate 
value of n within each interval, and the responses 
associated with the values of n within that in- 
terval are generated by the foregoing procedure. 
Note that higher accuracy of the reduced solu- 
tions can be obtained by marching backward as 
well as forward in the n-space with the reduced 
equations. 

3. The foregoing procedure can be directly applied to 
the solution of the governing nonlinear equations 
(eqs. (5)). This is accomplished by using a reduc- 
tion method, with the control parameter selected 
to be load, displacement, or arc length in the so- 
lution space and the global approximation vectors 
selected to be the various-order derivatives of the 
response quantities with respect to the control pa- 
rameter (see ref. 14). The global approximation 
vectors are obtained by successive differentiation 
of the governing equations (eqs. (5)) with respect 
to the control parameter. The left-side matrix of 
those equations has the same form as that of equa- 
tions (9). If the global approximation vectors are 
evaluated at zero value for the control parameter, 
the matrix {K}^'^ on the left side of equations (9) 
vanishes and the equations uncouple in harmon- 
ics. The application of the foregoing procedure 
considerably reduces the computational effort in 
generating the global approximation vectors and 
greatly enhances the effectiveness of the reduction 
method. 

The computational effort can be further reduced 
by using the procedure outlined in reference 14 
to uncouple the equations associated with the 
symmetric and antisymmetric shell parameters 
(with respect to 6 — 0). The procedure is based 
on transferring the anisotropic (nonorthotropic) 


([A'] + n„\K}) {AZ} rl „ ={P} + X{h (21) 

- [A]{AZ}„„ (22) 

(w + M»|)£{azU + 

-2[£]|-{AZ}„„ (23) 
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terms (submatrix [K a ] in eqs. (8)) to the right 

sides of equations (12) and adding another level 

of PCG iterations to account for them. 

Results and Discussion 

Numerical studies were performed to assess the 
accuracy of the two-dimensional shell model tire 
and the effectiveness of the computational proce- 
dure described in the preceding section, for generat- 
ing the response associated with different harmonics. 
Herein, the application of the model and the compu- 
tational procedure to the Space Shuttle orbiter nose- 
gear tire are presented. The geometric and material 
characteristics of the tire are given in figures 2 and 3. 
The Space Shuttle orbiter nose-gear tire is a 32 x 8.8 
type VII bias-ply tire with a ply rating of 16. The 
tire carcass is constructed of 10 lamina of nylon and 
rubber with an additional reinforcing ply beneath the 
tire tread as shown in figure 2. The tire has a three- 
groove tread, but the model assumes a smooth tread 
instead. The rated load for the tire is 15 000 lb at an 
inflation pressure of 320 psi. 

The numerical studies were performed with three- 
field mixed finite-element models used for the dis- 
cretization of the tire in the meridional direction. 
Linear interpolation functions were used to approx- 
imate each of the stress resultants and strain com- 
ponents, and quadratic Lagrangian interpolation 
functions were used to approximate each of the gen- 
eralized displacements. The integrals in the govern- 
ing equations were evaluated with a two-point Gauss- 
Legendre numerical quadrature formula. Because of 
the symmetry of the shell meridian and loading, only 
half of the tire meridian was analyzed. The finite- 
element models used are shown in figure 4. 

Modeling of the Tire Geometry 

The Space Shuttle orbiter nose-gear tire was mod- 
eled as a two-dimensional laminated shell with vari- 
able thickness and variable stiffness. The outer sur- 
face of the tire was taken to be the reference surface 
of the shell model. A tire was cut into sections and 
used to obtain accurate measurements for the cross- 
sectional shape of the uninflated tire. A smoothed 
spline under tension was used to fit a curve through 
the measured coordinates of the cross-sectional pro- 
file in a least-squares sense (see refs. 18 and 19). Be- 
cause of symmetry, only half of the cross section was 
modeled. A smooth variation of the second deriva- 
tive, “f , was achieved by adjusting the standard de- 
viations of the measured profile at the data points. 
For a detailed description of spline smoothing tech- 
niques, see references 19 to 21. 


The spline function, with the tension factor set 
equal to 0.1 and slope continuity enforced at both 
ends of the curves, was used to generate additional 
points along the tire meridian. The interpolation 
procedure is outlined in appendix A and the result- 
ing geometric characteristics of the tire are presented 
in appendix B and figure 5. The thickness of the 
tire carcass at the nodal points of the finite-element 
model was computed along the normal vector to the 
tire reference (outer) surface by locating the points 
of intersection of the normal vectors with the inner 
surface of the tire carcass. To facilitate these compu- 
tations the tire inner surface was approximated by a 
set of third-degree polynomials. 

Evaluation of Stiffness Coefficients of the 

Two-Dimensional Shell Model 

The cord-rubber composite was treated as a lam- 
inated material. For the purpose of computing stiff- 
ness variations in the meridional direction, the tire 
model was divided into seven regions, as shown in fig- 
ure 2. Thicknesses of the individual carcass plies were 
measured at the interfaces between the regions and 
are given in table 3. A linear variation was assumed 
for the thickness within each region. The thickness 
of the tire tread and sidewall covering was computed 
by subtracting the sum of the individual ply thick- 
nesses from the total thickness of the carcass at each 
location. 

The material properties of the different plies were 
obtained with the mechanics of material approach, 
which has been widely applied to rigid composites. 
(See refs. 22 and 23.) The elastic constants of the 
tire constituents used in this study are presented in 
table 4. It was assumed that nylon cords of two 
different diameters were used in the construction of 
the tire: d = 0.022 in. for the bottom two plies and 
the tread reinforcement in region I, and d = 0.031 in. 
for all other plies. 

The cord end counts (epi) for individual plies at 
the region interfaces are given in table 5. A linear 
variation was assumed for epi within each region. 
The formulas for evaluating the composite elastic 
coefficients for each ply, from the properties of the 
ply constituents, are given in appendix C. 

The stress-strain relations of the two-dimensional 
shell were obtained by first transforming the stiff- 
nesses of each of the individual layers to the global 
shell coordinates (s and 9) and then integrating these 
coefficients through the thickness. The cord orienta- 
tions in the individual plies of each region are given in 
table 6. The following formula was used to determine 
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0*., the angle (in degrees) measured from the s-axis 
to the 0-axis, at the numerical quadrature points: 

0 = Max { (54.382 - 3.884£ - 148.96£ 2 )° , 33° } 

(24) 

where £ is the dimensionless coordinate along the tire 
meridian. 

The resulting shell constitutive relations are given 
in appendix C. The meridional variations of the 
stiffness coefficients of the shell model are shown in 
figure 6. 


Case of Inflation Pressure 

To assess the accuracy of the shell model of the 
tire, the deformations produced by uniform inflation 
pressure of p Q = 320 psi, acting normal to the in- 
ner surface, were calculated using the geometrically 
nonlinear shell theory. Twelve finite elements were 
used in modeling half the cross section (a total of 
384 strain parameters, 384 stress- result ant parame- 
ters, and 243 nonzero generalized displacements; see 
fig. 4(a)). The calculated values were compared with 
the experimental data obtained on the Space Shuttle 
orbiter nose-gear tire. (See fig. 2.) The results are 
summarized in figures 7 to 9. Close agreement be- 
tween the predicted deformations and experimental 
results is demonstrated in figure 7. Figures 8 and 
9 show the meridional variations of the generalized 
displacements, stress resultants, and strain energy 
densities. As shown in figure 9 for the case of in- 
flation pressure, the transverse-shear strain energy 
density is considerably smaller than the extensional 
and bending strain energy density. 


Case of Localized Loading 

To assess the effectiveness of the computational 
procedure, linear solutions were obtained for a lo- 
calized normal loading on the outer surface simulat- 
ing contact pressure. The normal loading (in pounds 
per square inch) is given by the following equations, 
which model experimental data obtained at Langley 
on the shuttle tire: 


10 


P = 


where 


p n cosn0 (- 0.2 < £ < 0 . 2 ) 

7T - f 

n= 1 

o (|£|>o.2)J 


2 P° . n 
Pn = sin n(3 

717T 


> ( 25 ) 


(26) 


and p 0 and / 3 are functions of £ as shown in figure 10. 


Because of the symmetry of the shell meridian 
and loading, only half the meridian is analyzed us- 
ing 37 elements (a total of 1184 stress- resultant pa- 
rameters, 1184 strain parameters, and 743 nonzero 
displacement degrees of freedom; see fig. 4(b)). The 
boundary conditions at the centerline are taken to 
be the symmetric or antisymmetric conditions. Typ- 
ical results are presented in figures 11 and 12 and in 
tables 7 and 8. 

The foregoing procedure was applied to this prob- 
lem, and 10 global approximation vectors were eval- 
uated at n 0 = 5 and used to generate the tire re- 
sponse for n — 1 to 10. Accuracy of the generalized 
displacements obtained by the procedure with 8, 10, 
and 15 global approximation vectors is indicated in 
figures 11 and 12. Each generalized displacement in 
figures 11 and 12 is normalized by dividing by its 
maximum absolute value given in tables 7 and 8. 
Generalized displacements predicted by the foregoing 
procedure with 15 vectors are almost indistinguish- 
able from those predicted by the direct finite-element 
solution. 

Conclusions 

A computational procedure is presented for the 
geometrically nonlinear analysis of aircraft tires. 
The Space Shuttle orbiter nose-gear tire was mod- 
eled through use of a two-dimensional laminated 
anisotropic shell theory with the effects of variation 
in material and geometric parameters included. 

The governing discrete equations of the tire are 
obtained through application of the three-field Hu- 
Washizu mixed variational principle. The multilevel 
operator splitting is used to (1) uncouple the equa- 
tions associated with different harmonics, (2) identify 
the effects of different Fourier harmonics, and (3) de- 
lineate the effect of anisotropic (nonorthotropic) ma- 
terial properties. The nonlinear governing finite- 
element equations of the tire are solved with the 
Newton-Raphson iterative procedure. An efficient 
procedure is presented for the solution of the result- 
ing algebraic equations at each iteration, associated 
with different Fourier harmonics. The effectiveness 
of this procedure is demonstrated by means of a nu- 
merical example of the linear response of the Space 
Shuttle orbiter nose-gear tire subjected to inflation 
loading. The tire model is subjected to localized 
normal loading on the outer surface (simulating the 
contact pressure). 

Results of the present study suggest the follow- 
ing conclusions relative to the two-dimensional shell 
model used in simulating the response of the tire and 
to the proposed computational procedure for gen- 
erating the tire response associated with different 
Fourier harmonics: 
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1. A two-dimensional shell model with variable geo- 
metric and stiffness characteristics accurately pre- 
dicts the deformation of the tire when subjected 
to inflation pressure. 

2. Use of path derivatives (derivatives of the response 
with respect to the Fourier harmonic) as global 
approximation vectors leads to accurate solutions 
with a small number of vectors. Therefore, the 
time required to solve the reduced equations is 
relatively small and the total time required to 
generate the response for a range of 10 Fourier 
harmonics is little more than that required for a 
single Fourier harmonic. 

3. Global approximation vectors provide a direct 
measure of the sensitivity of the different response 
quantities to the circumferential wave (harmonic) 
number. Sensitivity of the global response can 
also be assessed with these vectors. 


4. The reduction method used in the proposed com- 
putational procedure exploits the best elements 
of the finite-element method and of the Bubnov- 
Galerkin technique, as follows: 

a. The finite-element method is used as a general 
approach for generating global approximation 
vectors. The full finite-element equations arc 
solved only for a single Fourier harmonic. 

b. The Bubnov-Galerkin technique is used as an 
efficient procedure for minimizing and distrib- 
uting the error throughout the structure. 

5. The reduction method extends the range of appli- 
cability of the Taylor series expansion by relaxing 
the requirement of using small changes in the cir- 
cumferential wave number. 

NASA Langley Research Center 
Hampton, VA 23665-5225 
January 22, 1990 
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Appendix A 


Summary of the Equations for a Spline 
Under Tension 

This appendix presents the interpolation proce- 
dure which uses a smoothed spline under tension to 
fit a curve through the measured coordinates of the 
cross-sectional profile in a least-squares sense. The 
cubic spline function requires that the second deriva- 
tive be piecewise linear and continuous. Thus, the 
second derivative with the effect of tension is speci- 
fied by 


x(s) - <7 2 x(s) = [x(.q) - g 2 y»]; _ * + [x(s i+1 ) - <7 2 y i+ il 

S l-f-l S l J 

(Al) 


• s i + 1 “ s i. 


where x = [x(s), 2 ( 5 )] is the position vector for points 
along the segment; y* = [x,z]i is the corresponding 
data at point i; a dot over a symbol refers to a 
derivative with respect to s. The chordal length 
(polygonal arc length) Sj is given by 


si — 0 


and 


(A2) 


Si — S Z _1 + [(xj — X-i-i) + (Zi — Z t -\) 


2i 1/2 


After solving equations (Al) for x(.s) and replacing 
x(sj) with y 2 , we obtain 


x(s) 


x(s,)/er 2 sinh — s) 


+ 


sinh<r(sj +1 — s 2 ) 
y i ~ x(s i)/a 2 ] (,s !+1 - s) 


s i 


+ 


x(s ?+ i)/cr 2 sinh a (s - Sj ) 
sinh + ] — s j) 

y?:+l - x(s J+ i)/cr 2 ] (s - s 

A ’? + l ~ s i 


(A3) 


From differentiating equations (A3) and equating 

right- and left-side derivatives at ,s t (for i = 2, 3 

/ — 1), we obtain a set of linear algebraic equations for 
x(s 2 ). With the assumption of a nonperiodic spline 
in which both slopes at s\ and s/ are provided, the 
tridiagonal differential equation is easily solved. 

Once the second derivatives at point i = 1 , 2, . . . , l 
are obtained, the first and second derivatives x(,s) 
and x(s) at the interpolated points arc evaluated 
by differentiating equations (A3) and they are used 
to compute the geometric parameters of the tire in 
appendix B. 

A normalized tension factor is used to eliminate 
a nonlinear behavior by setting (sec ref. 24) 


/ 


a 


°jsi - Si) 

l - 1 


(A4) 


In practice, if this factor is less than 0.001, the 
resulting curve is approximately a cubic spline, and 
if it is greater than 50, the curve is nearly piecewise 
linear. Note that s in equations (Al) to (A4) is not 
the actual arc length but the chordal length. Thus, 
the more data points one has from smoothing, the 
more accurate the arc length. 
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Appendix B 


Curvatures : 


Geometric Characteristics of the Tire 

This appendix presents the geometric charac- 
teristics derived from the application of the inter- 
polation procedure described in appendix A. From 
this spline interpolation, we obtain the following in- 
formation about the curve at each interpolated data 

point: af , a§» £ f > and The geometric 

parameters of the tire are then evaluated as follows: 


dx/ds > 

\/(dx/ds) 2 + ( dz/ds ) 2 

> 

dz/ds 

y/(dx/ds ) 2 4* {dz/ds) 2 


Normal vector. 


»=( sin *}= 
[ COS (j> J 


1 3 

K\ — — = -COS 0 

R\ 


( (fix 



K2= k = 


COS <f) 


Note that x = r, the normal distance from the axis 
to the reference surface, and ^ Also note that 

the normal vector is used to compute the thickness, 
the components of the inflation pressure acting on 
the inner surface of the tire, and the transformation 
matrix with respect to the global coordinate system, 
if necessary. 
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Appendix C 


Constitutive Relations for the Two-Dimensional Laminated Shell Model 

This appendix presents the formulas for evaluating the composite elastic coefficients for each ply (from the 
properties of the ply constituents) and the resulting shell constitutive relations. The stress-strain relations of 
the orthotropic and unidirectional layers are given by 


r o-l > 


2 

£ 

to 

1 


r ' 

cr 2 


C12 c 22 


£2 

^23 

► = 

C44 

< 

723 



C55 


713 

< T \2 > 


• C66 . 


• 712 > 


where the reduced stiffnesses Cjj are given by 


c 11 
C12 
c 22 


Ei 

1 - vi2V2\ 

E 1^21 
1 - 1 ^ 12^21 

E2 

1 - 1^12^21 


C44 = G23 


C55 = G13 


C66 = G12 


The elastic constants are computed by (see ref. 22) 


E c fc + E r ( 1 — f c ) 

Vcfc T v r ( 1 — fc) 

Er[E c ( 1 H- 2/,) + £ r (l - / r )] 

E c (\- f c ) + 2 E r {l+ 0.5fc) 

_ G r \G c + G r + (G c — G r )f r \ 

[G c + G r — (G r — G r )f c ] 

O.6G12 

V\2^2 

Ei 

where subscripts c and r represent the quantities of the nylon cord and the rubber, respectively, and f c is the 
volume fraction of the nylon cord: 

_ 7rd 2 (epi) 

/C 4 h k 

where d is cord diameter, h ^ is the layer thickness, and epi is the cord end count (in ends per inch). 


Ei = 
V 12 = 

E 2 = 

Gi2 = 
G23 = 
^21 - 
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The relationships between the stress resultants and strain measures of the tire are given by 


' N s 1 


A ii A\2 j 

Bn 

B 12 

© 

1 

1 




No 


a 22 @ j 

B 12 

B 22 

^26) 

1 

1 




N s () 


1 

^66 j 

© 

^26) 

Bm 

1 

1 



2 £ s0 

M s 

> = 


Dn 

Dn 

© 

1 

1 


< 

K s 

M 0 




D 22 

(S) 

1 



K 0 

MsO 






1 

1 



2k s q 

— 


Symmetric 



£>66 

1 



— 

Qs 






^55 

^54) 


2e S 3 

. Qo - 







^44 - 


. 2^03 , 


where and D v j (i, j = 1,2,6) are shell stiffness coefficients. The nonorthotropic (anisotropic) terms 

are circled and dots indicate zero terms. 
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Table 1. Symmetric and Antisymmetric Tire Parameters With Respect to 6 — 0 



Symmetric set 

Antisymmetric set 

Strain components 

{E}n 

£.s,n? ^,ni 

7i5 26^3^ j ^93 ,n 

£.s,7i! £0,ni ^£sO,m R'S.ni 

K().n i 2 K s 0,n> 2?s3,tn 2f^ 3 n 

Stress resultants 

{ H}n 

Nsjii N s Q n , AI S ji , 

-^.s0,n> Q#,™ 

Af s,ni AIq u1 N s q u , A / S) n? 
M 9, n ? A4 s Q n , Q$ji) Qo, 7i 

Generalized displacements 

Wn 

u ni Vui ^s,ni 4 > 9 m 

%, 0 t s,ni 00, n 


Table 2. Different Types of Coupling in Analysis of Tires With Semianalytic Finite Elements 


Response 

Material 

Governing finite-element equations 

Linear 

Isotropic or 
orthotropic 

Uncoupled in harmonics 

Symmetric and antisymmetric 
variables uncoupled 


Anisotropic 

Uncoupled in harmonics 



Symmetric and antisymmetric 
variables uncoupled 

Nonlinear 

Anisotropic 

Coupled in harmonics 



Symmetric and antisymmetric 
variables coupled 
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Table 3. Variation of Ply Thickness h k /h ( 
[h 0 = 0.7513 in.] 

h k jh 0 for region — 


and sidewall) 


0.0865- 

.0865 

.0865- 

.0865 

.0865- 

.0865 

.0666 

.0666 

.0666- 

.0666 

.0666- 

.0666 

.0666- 

.0666 

.0666- 

.0666 

.0666- 

.0666 

.0466- 

.0466 

.0466- 

.0466 

.0798- 

.0798 

0 


0.0865- 

.0658 

.0865- 

.0658 

.0666- 

.0506 

.0666- 

.0506 

.0666- 

.0506 

.0666- 

.0506 

.0666 

.0506 

.0666 

.0506 

.0466- 

.0354 

.0466- 

.0354 

.0798- 

.0798 

0 


0.0658- 

.0692 

.0658- 

.0692 

.0506- 

.0532 

.0506- 

.0532 

.0506- 

.0532 

.0506- 

.0532 

.0506- 

.0532 

.0506 

.0532 

.0354- 

.0373 

.0354- 

.0373 

.0798 

.0798 

0 


16 

= h- £ h k 

k = 2 

0.0692 

.0813 

.0692- 

.0813 

.0692 

.0813 

.0426 

.0500 

.0426 

.0500 

.0426 

.0500 

.0426 " 

.0500 
.0426 
.0500 
.0426 
.0500 
.0346- 
.0407 
.0346- 
.0407 
.0798 
.0798 

0 


0.0801- 

.0937 

.0801- 

.0937 

.0488- 

.0571 

.0488- 

.0571 

.0488 

.0571 

.0488 

.0571 

.0488- 

.0571 

“.0488- 

.0571 

.0488- 

.0571 

.0488- 

.0571 

.0375- 

.0439 

.0375- 

.0439 

.0798- 

.0798 

0 


0.0681- 

.1238 

.0681- 

.1238 

.0523- 

.0950 

.0523- 

.0950 

.0523 

.0950 

.0523- 

.0950 

.0523- 

.0950 

.0523 

.0950 

.0523- 

.0950 

.0523- 

.0950 

.0523- 

.0950 

.0523- 

.0950 

.0366- 

.0666 

.0366- 

.0666 

.0798- 

.1464 


0.0918- 
. 124C 
.0918- 
. 124C 
.0652- 
.088( 
.0652- 
.0881 
h .2662 
.359^ 
.0652 
.0881 
.0652- 
.088( 
.0652- 
.0881 
.0652- 
.0881 
*.2662 
.359 
.0652- 
.0881 
.0652- 
.088 
.0466- 
.062 
.0466- 
.062 
.1464- 
.159 


“Second layer of region I represents the layer which has the reinforcement (see fig. 2). 
6 This represents the thickness of the bead wires. 



Table 4. Values of Elastic Constants of Tire Constituents Used in Present Study 


Tire constituent 

Young's modulus, psi 

Shear modulus, psi 

Poisson’s ratio 

Rubber 

4.5 x 10 2 

1.51 x 10 2 

0.49 

Nylon cord 

3.5 x 10 5 

7.00 x 10 2 

.66 

Bead a 

2.9 x 10 7 

1.10 x 10 7 

.30 


Since the deformations are small in the bead area, it is reasonable to assume that the bead wires are 
isotropic. 


Table 5. Variation of Nylon Cord End Counts in Different Plies Along Meridian 


Ply number 
(top to bottom) 


Cord end count 

, ends per inch for region 



i 

ii 

III 

IV 

V 

VI 

VII 

1 tread and 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

sidewall 








2 

16-16 

18-14 

14-14 

14-14 

14-14 

14-14 

14 14 

3 

18 18 

18-14 

14-14 

14-14 

14 14 

14 14 

14-14 

4 

18-18 

21 20 

20- 18 

18-16 

16 16 

16-14 

14 14 

5 

23 21 

21-20 

20-18 

18 16 

16 16 

16-14 

14-14 

6 

23 21 

21 20 

20-18 

18 16 

16 16 

16-14 

Bead 

7 9 

23-21 

21-20 

20-18 

18-16 

16 16 

16-14 

14-14 

10 

23-21 

29-26 

26 25 

18 16 

16 16 

16-14 

14 14 

11 

30-29 

29 26 

26-25 

25-24 

16 16 

16-14 

Bead 

12 

30-29 

Rubber 

Rubber 

25 24 

24 22 

16 14 

14 14 

13 

Rubber 



Rubber 

24 22 

16-14 

14 14 

14 





Rubber 

22 22 

22 22 

15 






22-22 

22 22 

16 






Rubber 

Rubber 


18 





Table 6. Variation of Cord Orientation of Individual Plies, 0*-, Along Meridian Region 


0 k , deg. for region 


number 
o bottom) 

I 

ii 

hi 

IV 

V 

VI 

VII 

1 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

Rubber 

2 

-0-6 

0 

0 

0 

0 

0 

0 

3 

+0 

-0 

-0 

-0 

-0 

-0 

-0 

4 

-0 

0 

0 

-0 

-0 

-0 

0 

5 

0 

-0 

-0 

0 

0 

0 

-0 

6 

-0 

0 

0 

-0 

0 

-0 

Bead 

7 

0 

-0 

-0 

0 

-0 

0 

-0 

8 

-0 

0 

0 

-0 

0 

-0 

0 

9 

0 

-0 

-0 

0 

-0 

0 

-0 

10 

-0 

0 + 6 

0 + 6 

-0 

0 

-0 

0 

11 

0 + 6 

1 

1 

1 

1 

0 + 6 

-0 

0 

Bead 

12 

1 

Oil 

l 

Rubber 

Rubber 

1 

1 

0 + 6 

0 

0 

13 

Rubber 



Rubber 

-0-6 

-0 

-0 

14 





Rubber 

0 + 6 

04-6 

1 e: 






-0-6 

-0-6 

to 

16 






Rubber 

Rubber 



Table 7. Maximum Absolute Values of Normal Displacement Components w n 
[Tire shown in fig. 2; p 0 ,max = 445.3 psi; h a = 0.7513 in.; E To = 1160.3 psi] 


Fourier 
harmonic, n 

WnEj^/ (po.max^o) 

Fourier 
harmonic, n 

w n E r Y o / ( p 0 , max h a ) 

1 

18.08 

6 

1.095 

2 

11.41 

7 

.3122 

3 

7.110 

8 

1.134 

4 

4.841 

9 

1.424 

5 

2.906 

10 

1.313 


Table 8. Maximum Absolute Values of In-Plane Displacements and Rotation Components 
[Tire shown in fig. 2; p 0 , max = 445.3 psi; h () = 0.7513 in.; E r> = 1160.3 psi] 


Displacements 
and rotations 

Fourier harmonic 

n = 1 

n = 5 

n = 10 

U n Ej' o / (po, max h 0 ) 

7.584 

0.3266 

1.229 x 10- 1 

VnE^/ max h Q ) 

.04162 

.01139 

3.773 x 10 -3 

0s, tiEt o /Po, max 

3.602 

.4840 

1.304 x 10“ 1 

00.n^T fJ / Po, max 

.1805 

.01271 

6.982 x 10~ 3 


20 





21 


Generalized Stress resultants 
displacements 

Figure 1. Two-dimensional shell model of tire and sign convention for external loading, generalized displace- 
ments, and stress resultants. 
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Figure 3. Characteristics of Space Shuttle orbiter nose-gear 




Figure 4. Finite-element models used i 
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(a) Stiffness coefficients associated with uncoupled (orthotropic) response. 

Figure 6. Meridional variation of stiffness coefficients of two-dimensional shell model of Space Shuttle 
nose-gear tire. Ej^ o = 1160.3 psi. 
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Fig “ r ' 9 - Mer ; di »f 1 °f*™n energy densities produced by inflation pressure using two-dimensional 

shell model of Space Shuttle orbiter nose-gear tire. 
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Figure 11. Accuracy of normal displacements w n obtained with proposed procedure using two-dimensional 
shell model of Space Shuttle orbiter nose- gear tire. 
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Figure 12. Accuracy of in-plane displacements and rotation components obtained with proposed procedure 
using two-dimensional shell model of Space Shuttle orbiter nose-gear tire. 
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